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Abstract. Strain-induced deformations in graphene are predicted to give rise to large 
psendomagnetic fields. We examine theoretically the case of gas-inflated bubbles to determine 
whether signatures of such fields are present in the local density of states. Sharp-edged bubbles 
are found to induce Friedel-type oscillations which can envelope pseudo-Landau level features 
in certain regions of the bubble. However, bubbles which minimise interference effects are also 
unsuitable for pseudo-Landau level formation due to more spatially varying field proHles. 


1. Introduction 

Strain engineering has been proposed as a method to manipulate the electronic, optical and 
magnetic properties of graphene [T -10 .It is based on the close relation between the structural and 
electronic properties of graphene. An inhomogeneous strain field can introduce pseudomagnetic 
fields (PMFs), where the altered tight binding hoppings mimic the role of a gauge field 

in the low energy effective Dirac model of graphene 11,12 . Guinea et al demonstrated that 
nearly homogeneous PMFs can be generated by applying triaxial strain. One of the most striking 
consequences of homogeneous PMFs is the appearance of a Landau-like quantization. 
Scanning tunnelling spectroscopy on bubble-like deformations has observed pseudo-Landau levels 
corresponding to PMFs stronger than 300 T [^[^. 

Deformations can be induced in graphene samples by different techniques like pressurizing 
suspended graphene [5, 13 or by exploiting the thermal expansion coefficients of different 
substrates [^. As a result, introducing nonuniform strain distributions at the nanoscale is 
a promising route towards strain engineering. The standard theoretical approach to treat 
strain effects employs continuum mechanics to obtain the strain field. The strain field 
can then be coupled to an effective Dirac model of graphene to study the generation of 
PMFs in various geometries. In most studies, only the PMF distribution is considered as 
opposed to experimentally observable quantities like local density of states (LDOS). This study 
calculates the LDOS of such systems without applying periodicity, which can introduce spurious 
interactions between neighboring bubbles. 


2. Model 

2.1. Patched Green’s function approach 

The patched Green’s function approach, developed in Ref [^, treats device ‘patches’ embedded 
within an extended two dimensional system described by a tight-binding Hamiltonian. This 
approach allows us to insert a single bubble into an otherwise pristine infinite graphene sheet. 










Figure 1 . The patched GF method describes the extended graphene sheet away from the bubble device 
region (a) with a self energy term (b). We consider membrane and non-linear plate type bubbles 
with radially dependent (c) height and (d)strain profiles. 


and avoids issues such as interferences between a bubble and its periodic images or system edges. 
The extended part of the system is treated through a self-energy Tib entering the device area. 

Gd,d = {E1-Hd,d-^b)~\ ( 1 ) 


The self-energy is written in terms of Green’s functions (GFs) of an infinite, pristine sheet so that 
it can be calculated using methods taking advantage of periodicity and analytic integrability 15 


The Hamiltonian for the device region can be tridiagonalized allowing the GF of the device 
region, Gd,d, to be treated using an adaptive recursive method. A dual recursive sweep allows 
for efficient calculation of local properties everywhere in the device region surrounding a bubble, 
enabling us to investigate spatial variations of the LDOS. 


2.2. Strain model 

In this work we consider two possibilities for the shape of a gas inflated bubble, namely 
the membrane and non-linear plate models. The membrane model is suitable for very 
large bubbles where bending stiffness can be neglected, whereas the non-linear plate model 
is more appropriate for including bending effects near the edges of smaller bubbles. Membrane 
bubbles therefore have very sharp edges, whereas the edges are smoother in the non-linear plate 
bubbles. While these continuum models have been found to agree well with experimental shape 
profiles, more accurate modeling of bubble shapes and strain distributions can be achieved using 
molecular dynamics simulations. (5 ,17 

Using the deformation field u = {ur, ug, z) the position of the atom i initially at becomes 
Ri = R^ + u. The new bond lengths are afterwards determined as dij = \Ri — Rj\. Hoppings 

are modified according to Uj = where to « 2.7 eV is the pristine coupling, 

oo = I. 42 A is the carbon-carbon distance and /3 ^ 5log(t)/(9log(a)|a=ai 3 ~ strain 

tensor created by the strain field, eij = -|- d^uj -\- d.^zdjz), gives rise to a gauge field |ll]. 


A = - 


hjd 

2eao 


c.T.T e 


yy 


2e 


xy 


( 2 ) 


with the resulting PMF given by Bg = V x A. Table shows the height profile, in-plane 
displacements and PMFs for the two types of bubbles considered here. 






















Model type 

z{r, 9) 

II 

Bs{r,9) 

Membrane 

1 

1—1 

O 



Non-linear plate 


/ r{R - r)(ci -1- C 2 r)\ 

V 0 ) 

32/i2rV r2\l 

(l ^,) sm(3») 


Table 1. Height profile z{r,9), in-plane displacements u{r,0) and PMF distributions Bs{r,6) 
for the membrane and non-linear plate bubbles, where R and /iq are the bubble radius and 
height respectively, and uq = 1.136/ig/ii, ci = l-SOS/ig/iZ^ and ci = —l.dSlh'^/R'^. 


3. Results 

New hopping parameters, calculated from atomic coordinates generated by the continuum model 
displacements, give a tight-binding description of the bubble region which can be used within the 
patched GF approach. From this, the LDOS at every site in the bubble, as well as the average 
DOS, can be quickly calculated. Fig. shows the averaged DOS for membrane (blue) and non¬ 
linear plate (red) model bubbles. For the membrane model, we have previously distinguished 
between two different type of oscillations [^. A series of sharp peaks, such as those highlighted 
by the blue circle and triangle, were found to have an energy dependence, En ~ ^/n, consistent 
with Landau-like levels arising due to the PMF. In addition, the periodic oscillations visible at 
higher energies are identified as Friedel-type oscillations arising due to scattering of electrons 
induced by the sharp edges of the membrane-model bubble. 

Before discussing the interplay between different oscillation types, we note that both features 
vary independently with position throughout the bubble region. The position dependence of the 
pseudo-Landau level features arises due to the non-uniform PMF distribution within the bubble, 
which is plotted in Fig. [^. This takes maximum amplitudes along the armchair directions which 
occur every 60°, but is only three-fold symmetric due to a sign change between two consecutive 
amplitude maxima. Unlike real magnetic fields, PMFs conserve time-reveral symmetry by taking 
opposite signs in the K and K' valleys of graphene. One manifestation of this is a strong 
sublattice polarisation (l4, 19 20 , which is clearly visible for the sublattice-split LDOS maps 


shown in Fig. [^-c, where the circle LDOS peak from Fig. is localized in different regions for 
different sublattices. Comparison to panel confirm that these correspond to a change in the 
sign of the PMF. We note that this, the first “pseudomagnetic peak”, is localized along armchair 
directions where the PMFs are largest and reasonably constant. The position dependence of the 
Friedel oscillations meanwhile emerges from interference between electrons scattered at different 
sides of the bubble. The interplay of both oscillation types in membrane bubbles leads to the 
Friedel type acting as an envelope and quenching the LDOS signature of pseudo-Landau peaks 
in certain regions of the bubble, as is clear for the dark spots in the LDOS map for the higher 
energy in Fig. Hi- STS measurements taken at such a spot would completely omit this peak due 
to the enveloping effects of the edge-induced Friedel oscillations. 

When searching for signatures of PMFs in gas-inflated bubbles, it may thus be worth 
considering bubbles with a softer edge profile, such as the non-linear plate model, which should 
give riser to weaker Friedel oscillations. This is clear from the averaged LDOS curve in Fig. 
[^, where the higher energy oscillations are considerably suppressed compared to the membrane 
case. However, in this case there is also an absence of sharp Landau-level-like peaks following 
a ^/n distribution, with the possible exception of the peak denoted by the red circle. This lack 


















Figure 2. a) shows the averaged DOS with each bubble model with R = lOnm and ho = lOnm. 
Important peaks in each are highlighted by symbols. b,c) The LDOS for the membrane 
model at the circle energy is mapped for the A and B sublattices separately. d,e) show the 
PMF distributions for each bubble type and f-k) show A sublattice LDOS maps for the peaks 
highlighted in a). The scale bar in all LDOS maps is 5nm. 


of pseudo-Landau features is consistent with the PMF distribution for this bubble type, plotted 
in Fig. [^. We note that this bubble, with less sharp edges, also has a radial fluctuation in the 
sign and strength of the PMF. The center of the bubble has a field distribution similar to that 
of the membrane case, and the central region of the LDOS map in Fig[^ resembles that of the 
corresponding membrane model peak (Fig[^). We note also that the Friedel features for the 
higher (triangle) energy in Fig. (which does not fit the ^/n distribution) are more blurred than 
for the membrane case, as expected for scattering from a less-sharp bubble edge. Thus it seems 
that bubble shapes which reduce Friedel oscillations also effectively remove pseudomagnetic 
Landau effects due to the less uniform PMFs induced by their strain profiles. 

Finally, we note that the square symbol energy peak at low energies in both bubble types is 
a state localized near the bubble edge and plotted in Figs. E and[^. It is not directly related to 
pseudomagnetic effects, but emerges due to the interface between the pristine graphene region 
outside the bubble and the strained, perturbed region within. The presence of localized states at 
this boundary acts somewhat like a potential, and induces the scattering which lies behind the 
Friedel oscillations in these bubbles. We note that these states in the non-linear plate bubble are 
far less localized than their membrane bubble counterparts, due to an edge which is no longer 
as sharp. This in turn leads to the smoothening and averaging out of the Friedel oscillations 
that we observed earlier for the non-linear plate bubbles. 

4. Conclusions 

We studied theoretically the local and averaged densities of states in gas-inflated graphene 
bubbles embedded in infinite graphene sheets by making use of the patched Green’s function 
approach. We determined that pseudo-Landau level features in sharp-edged bubbles may be 
hidden by interference effects due to electron scattering at the bubble edges. Softer-edged 
bubbles were found to display weaker interference effects, however their shape profiles also 
resulted in pseudomagnetic field distributions unsuitable for pseudo-Landau level formation. 
Our results suggest that it will be difficult to obtain reliable Landau level features in such gas 
inflated systems, unlike bubbles formed on substrates which often display the triaxial-type strain 
which is predicted to give a more appropriate pseudomagnetic field for Landau level formation. 
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